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We study the effective parameters of QED near the decoupling temperature and show that the 
QED perturbation theory works perfectly fine at temperatures, below the decoupling temperature. 
Temperature dependent selfmass of electron, at T = m, gives two different values when approached 
to the same overlapping point. It ia shown that at T = m, change in thermal contribution of the 
electron selfmass is 1/3 of the low temperature value and 1/2 of the high temperature value. The 
difference of behavior measures the electron background contributions at T = m. These electrons 
are emitted through beta decay. This rise in mass affects the QED parameters and change the 
electromagnetic properties of the medium with temperature also. However, these contributions are 
ignorable near the decoupling temperature. 
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Renormalization techniques of perturbation theory are used to calculate temperature dependence of renormalization 
constants of QED (quantum electrodynamics) at finite temperature [1-19]. The values of electron mass, charge and 
' wavefunction, at a given temperature, represent the effective parameters of QED at those temperatures. The magnetic 
moment of electron [9,18], dynamically generated mass of photon [4,9,15,19] and QED coupling constants are estimated 
' as functions of temperature. Moreover, thermal contributions to the electric permittivity, magnetic permeability and 
, dielectric constant of a medium can be obtained from the photon selfmass. Some of the important parameters of 
£ — ■ QED plasma such as Debye shielding length, plasma frequency and phase transitions can also be determined from 
C*~) ' the properties of the medium itself. 

In this paper, we quantitatively analyze the existing results of temperature dependent renomalization constants. 
The renormalization scheme of QED in real-time formalism, is used to calculate the electron mass, wavefunction 
and charge of electron. It is now well-known that the existing first order thermal corrections to the renormalization 
constants give the quadratic dependence of QED parameters on temperature T, expressed in units of electron mass m. 
The existing scheme of calculations is very useful and it works perfectly fine below the decoupling temperature, i.e.; 
2MeV. It can be explicitly checked that in the existing scheme of calculation, QED theory remains renormalizeable 
at T < 5MeV. However, the first order corrections will exceed the original values of QED parameters at much 
higher temperatures and hard thermal loops will not appear before that. Then we have to look for a new scheme of 
calculations for temperatures much higher than the decoupling temperatures. There are already developed methods 
[20,21], which could be used for this purpose. However, below the neutrino decoupling temperature, the real part of 
the propagators serve the purpose. 

Existing analytical results are based on the methods of perturbation theory in vacuum and are extremely useful 
to explain the behavior of QED up to the decoupling temperature. Renormalization methods of QED at finite 
temperature ensures a divergence free QED in a thermal medium below the decoupling temperature. Without going 
in to the calculational details, we give a brief overview of the existing results in the real-time formalism. It is possible 
to separate out the temperature dependent contributions from the vacuum contribution as the statistical distribution 
functions contribute additional statistical terms, both to fermion and boson propagators, in the form of Fermi-Dirac 
distribution and Bose-Einstein distribution functions, respectively. The Feynman rules of vacuum theory are used 
with the statistically corrected propagators given as 
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for bosons, and 



'Electronic address: masood@uhcl.edu 



2 



Sf(p) = 
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(2) 



for fermions. 

The electron mass, wavefunction and charge are then calculated in a statistical medium using Feynman rules of QED, 
with the modified propagators given in Eqs.(l) and (2). Since the temperature corrections are additive corrections 
in the propagator and appear as additive terms in the matrix element, thermal radiative corrections can be studied 
independent of vacuum corrections at the one loop level. We restrict ourselves to the one loop contributions only as it 
can be easily shown that the higher order contributions [13-19] are smaller than the first order contributions at these 
temperatures and ensure the validity of the renormalization scheme. 

II. SELFMASS OF ELECTRON 

The renormalized mass of electrons to# can be represented as a physical mass m p h ys oi electron and is defined in a 
hot and dense medium as, 



rriR = niphys = m + 5m(T = 0) + Sm(T). 



(3) 



where to is the electron mass at zero temperature and 5m(T = 0) represents the radiative corrections from vacuum 
and Sm(T) are the contributions from the statistical background at nonzero temperature T. The physical mass can 
get radiative corrections at different orders of a and can be written as: 



m phys = to + 5m {l) + <5to (2) + .... 



(4) 



where 5m^ and Sm^ are the shifts in the electron mass in the first and second order in a, respectively. The 



physical mass is deduced by locating the pole of the fcrmion propagator 



i(y+m) 



in thermal background. For this 



p 2 — m 2 +i£ 

purpose, we sum over all the same order diagrams. Renormalization is established by demonstrating the order-by- 
order cancellation of singularities. All the finite terms from the same order in a are combined together to evaluate 
the same order contribution to the physical mass given in cq.(4). The physical mass in thermal background, up to 
order a 2 [13-19], is calculated using the renormalization techniques of QED. Writing the selfmass of electron [1,6] as: 



E(p) - A{p)E lo - B(p)p.j - C(p), 



(5) 



where A(p), B(p), and C(p) are the relevant coefficients that are functions of electron momentum only. Taking the 
inverse of the propagator with momentum and mass terms are separated as: 



S _1 (p) = (1 - A)E 1 ° - (1 - B)p. 1 - (m - C). 



(6) 



The temperature-dependent radiative corrections to the electron mass up to the first order in a, are obtained from 
the temperature dependent propagators as 



m 2 phys = E 2 ~ |p| 2 = m 2 



6cn 

1 - — b(mp) 

TV 



+ 



4a 



mT a(m/3) + -airT 2 ~c(m/3) 

3 ir z 



(7) 



giving 
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where — is the relative shift in electron mass due to finite temperature which was determined in Ref. [7] with 



o(m/3) =ln(l + e- m/3 ) 



(9) 
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b(m0) = J^(-l) n Ei(-nm0), (10) 



n=l 



c(m/3) = (11) 

n=l 

The validity of Eq.(4) can be ensured for T < 2MeV, the decoupling temperature, as ^ is always smaller than 
unity within this limit. This scheme of calculations will not work for higher temperatures, as the summation over all 
orders of perturbative correction may exceed the original values of QED parameters, after 5MeV. We will have to 
develop a new method of calculations for the higher temperature limit. At low temperature T < to, the functions 
a(m(3), b(m/3), and c(m/3) fall off in powers of e _m/3 in comparison with (^) and can be neglected in the low 
temperature limit giving, 

Sm rrr , . anT 2 ... 
^ T< ^ = ^?- (12) 
In the high-temperature limit, a(m(3) and b(m(3) are vanishingly small whereas c(m(3) — ► — 7r 2 /12, yielding 

Sm frri . airT 2 

m (T>m) = ^- (13) 

Eqs.(12) and (13) give ^ = 7.647 x 1CT 3 ^ for low temperature and ^ = 1.147 x 10" 2 ^ for high temperature 
showing that the rate of change of Sm is larger at T > to as compared to T < to. Subtracting eq.(12) from (13), the 
change in ^ between low and high temperature ranges can be written as 

TO DTO TO 

showing that the A(^) = ±3.8 x 10~ 3 at T = to. It can be easily checked that the low temperature behavior will give 
a 50% decrease in selfmass as compared to high T value. Whereas the high T behavior will give 33% more selfmass as 
compared to the low T value. This difference is due to the fact that at low temperature, only hot boson contribution 
is calculated using Eq.(12), whereas Eq.(13) includes the fermion background contribution also. Eq.(14) is a measure 
of fermion background contribution that is ignoreable at T < to but cannot be ignored at T > to. This difference 
keeps on increasing with temperature also. 

Temperature dependence of QED parameters is a little more complicated and significant because of the change in 
matter composition, during nucleosynthesis. Therefore Eqs.(7) and (8) are required for T ~ to region and help to 
compute the change in thermal behavior of QED parameters due to the change in matter composition, carried out by 
beta decay and other processes, at that time. This difference can be clearly seen in Figure(l). We plot of eqs. (12) and 
(13), corresponding to low temperature and high temperature and show that both plots start to give a disconnected 
region near T <~ to, i.e; the nucleosynthesis temperature. Slope of both graphs is also different, indicating that T ~ m 
induces a change in thermal properties for a heating and a cooling system. 

The disconnected region in Figure (1) shows that the ^ behavior is needed to be modified to find the missing link 
between the low temperature and high temperature behavior. Eq.(8) provides the information about the disonnected 
region and eqs. (12) and (13) can be derived from eq.(8). Figure (1) also shows that the thermal corrections to the 
electron mass are not significant below T <~ 0.05MeV. 

The modification in the electron mass behavior in the range T <~ to, is estimated by Eq. (12). It is also clear from 
Eq. (13) that after 5MeV, temperature dependence correction term can easily reach (— > 1), even at the one-loop 
level. Higher order corrections [13-16] will grow up rapidly at high temperatures as we can approximate it as: 



m>piiy S — Tn 



Sm 1.<5to, 9 
m + 2 { m ) + 



Sm 

TOexp( ) (15) 

TO 



for the temperature T in the units of the elctron mass m as O.Ito <T> IOto. This almost exponential growth in 
the physical mass is very important here. 

Figure (1) shows that a change in the QED behavior occurs around T ~ to and is clearly related to nucleosynthesis. 
Right after decoupling, beta decay and other processes involving the electron mass, change the composition of matter 
and electron picks up thermal mass from hot fermion loop. At T > 5 MeV, the renormalization scheme of perturbative 
QED may not be valid as beta decay contribute through weak interactions. 
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FIG. 1: Low temperature values of selfmass of electron (solid line) and high temperature values (broken line) are plotted as a 
function of temperature T. Around T=m, both of the graphs approach to different values. 

III. WAVEFUNCTION RENORMALIZATION: 

The electron wavefunction in QED is related to the selfmass of electron through Ward identity. The factor (1 — A) 
is required for renormalization, because then the propagator can also be renormalized by replacing 



fa — m + ie fa — m + ie 

Thus, for Lorentz invariant self-energy, the wavefunction renormalization constant can equivalently be expressed as 

Z^ 1 = 1 - A 



J 2 

z 2 = l -^j- ( 16 ) 

The fermion wavefunction renormalization in the finite temperature field theory can be obtained in a similar way as 
discussed in vacuum theories. However, the Lorentz invariance in the finite temperature theory is imposed by setting 
A = B in Eq.(5). Thus, using Eqs. (16) and (5), one obtains [7] 



Z?(mfi) = Z-\T = 0)-- [°° ^n B {k) - ^6(m/3) 

7T J k 7T 



irvE 2 1 — V 6 
giving the low temperature values as 



aT 2 1 + v it 2 

In {— +m/3a{ml3) - c(m(3)}. (17) 



i i , , 2a f dk , x anT 2 1 1 - v 

and high temperature value as 



6E 2 v 1 + v 



(18) 



^ = = 0)-=J fn B (k) + ^-In^ (19) 

For small values of v, the low and high temperature values can be determined from eqs. (18) and (19) as 

„ i „ i n 2a f dk ,, , anT 2 

Z^ 1 = Z^(T = 0) - — J T n B (k) - (20) 



2a f dk anT 2 1 . 1- v 



■5 



for low temperature, and 



V = V(T = 0) -£/fn fl (*) (21) 

for high temperature. 

The finite part of eqs. (20) and (21) is equal to ^ at the lowest value of energy, i.e.; E = to in the relevant 
temperature range. These terms are suppressed at large value of electron energy E as they are suppressed by a 
factor Jy. However, the calculated value at that temperature is significantly different. The difference in the thermal 
contribution can easily be found to be about 50% of the low temperature value and around 33.3% of the high 
temperature value. This difference can be mentioned as 



A ^ r 1 )«-^ F = 3.8x10 (22) 

and is similar to the selfmass corrections. However, a difference of sign is noticeable, showing the decrease in the 
renormalization constant and not the increase. However, the finite part of the wavefunction renormalization constant 
can be obtained by finding a ratio of temperature with the Lorentz energy E. The minimum value of this energy is 
equal to mass. Following eq.(15), the higher order contributions to the wavefunction can then be approximated as 



i 1 , 2a f dk , . 

Z^ 1 = Z^(T = 0) - — / — n B (k) - exp( 



OtTtT 

~2E? 



(23) 



As shown in Figure 2, the wavefunction renormalization contributions are negative. It shows that the low temperature 
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FIG. 2: Low temperature value (broken line) and high temperature value (solid line) of the electron wavefunction renormal- 
ization constant are plotted as a function of T. 

(T < to) contribution is simply ignorable as compared to high temperature (T > m) contribution everywhere below 
the decoupling temperature. The finite term can be ignorable at T — to, and even at the temperatures higher than 
nucleosynthesis as T < E at those temperatures because E is always greater than m. For large E, thermal contributions 
are even smaller and more ignorable. So the two interesting physical limits give smaller thermal contribution in electron 
wavefunction as the relevant temperature limits can be defined as to < T < E and T < to < E, which ensures the 
renormalizability of QED at comparatively higher temperature as compared to selfmass. 

We do not discuss the infrared singularity term as it has already been studied in literature. We are interested 
to compare the finite contributions only to see how they change in different regions of temperatures can justify the 
validity of renormaliztion scheme of vacuum theories. 
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IV. PHOTON SELFMASS AND QED COUPLING CONSTANT 

Selfmass of photon and the electron charge also behave differently for a cooling and a heating system at high 
temperature. The electron charge is not expected to be changed due to the presence of neutral photons. That is 
the reason that the electron charge and the coupling constant does not show significant temperature dependence for 
T < m. However, they have significant thermal contributions at high temperatures (T > m). Difference of behavior of 
a cooling and a heating system starts to be significant near T ~ m from the low temperature side and near decoupling 
temperature from the high temperature side. This difference in the coupling constant occurs due to the dynamically 
generated mass of photon that couples with the hot electrons in thermal medium. The emission of electrons looks 
more natural due to the beta decay processes during nucleosynthesis. 

Calculations of the vacuum polarization tensor [7] show that the real part of the longitudinal and transverse 
polarization components of the polarization tensors can be evaluated, in the limit to — > 0, as: 



Real (M) = ^ {t 2 + , Jien* (*,o) = | f* 9 ^ + } • 

giving the interaction potential, in the rest frame of the charged particles as 



(24) 



V(k) 



2 c 
e B<V0 



+ 



k 2_ 4|a| T 2 + ^ ln m| fc 



9}±v y>fj,y>i/ 



(25) 



V(k) = efjc^oA^^o, 

en is the renormalized charge at T = 0. V(k) can be expanded, at low temperature as: 



(26) 



k 2 



iizaT 2 



k 2 



(27) 



The constant in the longitudinal propagator is the plasma screening mass, therefore, whole of the outside factor 
corresponds to the charge renormalization and in turn to the coupling constant. We may then write the coupling 
constant at low temperature as [1-3]: 



a (T) = a (T = 0) 1 



2a , T 
— In — 

Sir m 



a(T = Q) 1 + 1.55 x l(T 3 m- 



(28) 



The factor 1.55 x 10~ 3 In ^ is a slowly varying function of temperature and does not give any significant contribution 
near the decoupling temperature and remains insignificant for a large range of temperature. 

The temperature dependent factor in the longitudinal propagator 47r ° T is the plasma screening frequencies or 
selfmass of photon that contribute to the QED coupling constant at finite temperature. 

For generalized temperatures, the charge renormalization constant Z 3 can be written as [9] 



2e 2 ,c(m/3) ma(m(3) 1 2 w 2 ,,, „ N1 



P 2 



13 



(29) 



Also, the electric permittivity is 

4e 2 / .2 



e{K) ~ 1 + 



ir 2 K 2 



1 



CO' 



1 



uj uj + k 
— In 

2k uj — k 



and the magnetic permeability is 

1 2e 2 



1 



li(K) n 2 k 2 K 2 ' 
' c(m/3) ma(mj3) 



2 



c{mj3) ma(m/3) 



P 2 



k 2 



6m 



/3 



1 fc2 

UJ 2 



2 uj 2 



12% 



l/ 2 2 llfc 2 + 37w 2 \ , 
- f 2m 2 - uj 2 + ) 6(m/3)}, 



(30) 



1 — —r- I — - In 



to , w + fc 



fc 2 7 2fc w 
109fc 2 



72 



6(m/3)]. 



(31) 
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Electron Mass 


Charge 


Wavefunction for E=r 




LowT 


High T 


HighT 


LowT 


High T 


0.02 


3E-06 


4.59E-06 


4.9E-07 


-3.05E-06 


-4.6E-05 


:.C2 


o.oos 


0.011933 


0.00127 


-0.007939 


■c.c::9i 


2.02 


0.0312 


0.0+6302 


0.00495 


■C.C31136 


-0.0457 


3.02 


0.0697 


0.104611 


ecu: 


-0.069595 


-C.1CX39 


-.02 


0.1236 


0.13536 


0.01966 


-0.123315 


-0.13497 


5.02 


0.1927 


0.239049 


0.03066 


-0.192297 


-0.23345 


6.02 


0.2771 


0.415677 


C.C-^C9 


-0.27654 


3. 


7.02 


0.3763 


0.565246 


0.05995 


-0.376044 


-0.56407 


a. 02 


0.4919 


0.737755 


0.07325 


-0.49031 


-0.73621 


902 


0.6222 


0.933204 


0.09393 


-0.620337 


-0.93126 


1C.C2 


0.7673 


1.151593 


0.1221^ 


-0.766126 


-1.14919 


11.C2 


0.9237 


1.392921 


0.14774 


-0.926675 


-1.39001 



TABLE I: Renormalization constants of QED at finite temperature corresponding to the T < m and T > m expressions. 



In the limit T > m, the wavefunction renormalization constant can be written as: 



6m 2 



(32) 



giving the renormalized coupling constant as 



47T£o 6m 2 2h 6m 2 ' 

Eq. (32) gives Zs = 1 + 1.2 X 10 -3 ^ and leaves the perturbation series valid for at least T < 10m. 

It is clear from eqs. (29-32) that the coupling constant is basically changed from the hot fermion loop contributions. 
Hot bosons do not change the coupling constant and the vacuum fluctuations occur due to fermion loops at the first 
order in a. 

Temperature corrections to the coupling constant start to become noticeable during nucleosynthesis. However, the 
value of c(mf3) approaches a constant value = rw- at T ~ m and start to become significant at higher temperatures. 
After T > 5MeV, thermal contributions indicate that the coupling constant can grow larger than unity at high 
temperature indicating a problem for perturbative behavior of QED. We need to use non-perturbative methods to 
establish renormalization of QED at those temperatures. 



V. RESULTS AND DISCUSSION 



Quantitative study of renormalization constants at finite temperature [Table 1] shows that all the renormalization 
constants are definitely finite below 5MeV (i.e.;T < 10m). Table 1 indicates that the first order corrections to all the 
renormalization constants of QED are not always ignorable, and significantly grows around the decoupling temper- 
ature. Computation of thermal contributions of the renormalization constants, around the decoupling temperature, 
is not straightforward. In this range, the largest thermal contribution comes from the electron sclfmass. However, 
thermal contribution to the wavefunction renormalization constant is small at these temperatures, and high energies 
[E >> m). Eqs. (18) and (19) show that the maximum thermal contribution to electron wavefunction renormal- 
ization constant is equal to the selfmass of electron at E = m and that is the largest contribution. There is no 
low temperature contributions to electron charge, as well as the QED coupling constant, due to the absence of hot 
fermions in the medium. Hot fermion loop contributions are ignorable at low temperatures (T < m). However, at 
T > m the coupling constant starts to pick up thermal corrections, due to the dynamically generated mass of photon, 
which is generated through the fermion background only. 
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FIG. 3: Low temperature values of the renormalization constants of QED (solid line) and high temperature values (broken 
line) are plotted as a function of temperature T. Around low temperatures, both of the graphs approach to different values. 

A comparison of the statistical background contributions to different renormalization constants shows that, in 
this scheme of calculations, the selfmass corrections are the largest corrections at all temperatures. Wavefunction 
renormalization is significant at low energies but the coupling constant does not seem to get significant thermal 
corrections at decoupling temperature. This scheme of calculation is very helpful to compute QED parameters during 
nucleosynthesis in terms of a(m/3), b(m/3) and c(m/3). Figure 3 gives a plot of thermal contributions to electron mass 
), wavefunction renormalization constant (Z^ 1 + / jrifl(i) ) at low energy {E = \f T 2.rn) and electron charge 
(Z?, — 1 ), that can be derived for T << m and T » m ranges from the same equation. Since we are dealing with the 
exponential functions in this study, an order of magnitude difference is a safe limit for these approximations. It means 
T < 0.1m is T << m and T > 10m is T » m. The broken line is high temperature limit (T ^ m ) and the solid 
lines correspond to the low temperature (T < m ) limit of the corresponding parameters. Both limits are plotted for 
overlapping temperatures to show that the difference between the low temperature and the high temperature values 
is due to the contributions of hot fermion loops at high temperatures. Contributions of function c(mf3) vanishes 

at low T and it sums up to (-yj) f° r l ar S e T values. This function actually bring the fermion loop contributions 
at high T when more fermions are generated during the nucleosynthesis and their presence is not ignorable in the 
medium. Separations between two curves in the above graphs of Figure 3, measure the fermion contributions at high 
temperatures. However, around the decoupling temperature, this separation indicate complicated processes which 
cause the emission of electrons in beta decay and even absorption, during nucleosynthesis. These electrons acquire 
thermal equilibrium with the medium. However, a(m(3), b(mf3) and c(m/3) functions are needed to study the behavior 
of these parameters near T ~ m and explain the disconnected region of Figure (1). 

Plot of these renormalization constants, show that the temperature corrections are to small to differentiate between 
the boson and fermion contributions, for T sufficiently smaller than m. Significant thermal behavior starts near 
T ~ m, as indicated in Figures 3. 

Cooling universe of the standard big bang model behaves differently after the neutrino decoupling. Nucleosynthesis 
starts right after the neutrino decoupling and the helium synthesis takes place when the temperature of the universe is 
cooled down to the temperature of electron mass. This is actually a period, where the finite temperature corrections 
to QED parameters are significant and complicated enough to evaluate it numerically. However, the temperature 
dependent QED parameters are helpful to describe the observations of WMAP (Wilkinson Microwave Anisotropy 
Probe) data [22-23]. After the nucleosynthesis is complete, all the renormalization constants depend quadratically 
on temperature, though it may not always be significant. At high temperatures, QED coupling plays its role in 
modifying QED parameters for nucleosynthesis. With the help of these effective parameters of QED, the abundance 
of helium in the early universe can be estimated [12] precisely at a given temperature. The temperature dependent 
QED corrections to the nucleosynthesis parameters improve the results of standard big bang model of cosmology and 
testing of the standard model with WMAP becomes more reasonable. The same techniques can even be used to 
calculate perturbative effects in QCD [24] and electroweak processes [25] at low temperatures. 
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